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that  an  attenuating  relaxation  coefficient  is  recommended.  However,  when  wind  observations 
at  the  outflow  level  are  added  into  the  assimilation,  forecasts  improve  substantially.  Best 
forecasts  are  achieved  when  observations  over  the  entire  lower  troposphere  are  assimilated. 

Additional  experiments  indicate  the  errors  in  the  satellite  observation  contaminate  the 
forecast.  But  the  assimilation  of  inflow  and  outflow  winds  still  improve  the  intensity  forecast  if  the 
satellite  observation  errors  are  less  than  or  about  the  same  magnitude  of  those  in  the  initial  wind 
field. 
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A  SIMULATION  STUDY  ON  THE  IMPACT  OF  SATELLITE-SENSED  WINDS 
ON  TROPICAL  CYCLONE  FORECAST 

1.  Introduction 

The  improvement  of  forecast  skill  on  tropical  cyclones  evident  in 
the  1960's  has  not  been  continued  in  the  1970's  in  spite  of  improved  tech¬ 
nology  and  continuing  effort.  The  lack  of  improvement  has  been  attri¬ 
buted  to  the  imperfect  knowledge  of  the  initial  fields  for  objective  models. 
Elsberry  (1977)  attributed  the  poor  performance  of  his  prediction  model 
to  the  deficient..'  of  the  initial  wind  data.  For  the  1976  Atlantic  tropical 
cyclone  season,  Hovermale  and  Livezey  (1977)  showed  the  errors  for  the 
36  h  and  48  h  forecasts  increased  by  approximately  a  factor  of  three  for 
storms  over  data-void  ocean  regions  as  compared  to  storms  near  coastal 
stations.  In  addition,  the  theory  of  geostrophic  adjustment  requires 
that  the  mass  field  adjusts  to  the  momentum  field  for  low  latitudes  and 
tropical  cyclone  scales  of  motion  (Monin  and  Obukhov,  1959;  Washington, 
1964).  It  is  beyond  doubt  that  wind  observations  are  essential  to  tropical 
cyclone  forecasts. 

Although  initial  wind  analyses  of  tropical  cyclones  have  been 
-improved  from  reconnaissance  aircraft  flights,  the  quantity  and  especially 
the  aerial  coverage  of  wind  data  so  obtained  are  inadequate  for  numerical 
no  el  initializations.  Remote  measurements  from  geostationary  and  orbiting 
satellites  will  be  relied  upon  as  important  data  sources.  Rodgers  et  al 
(1979)  explained  technioues  to  derive  low-level  and  outflow  level  winds 
for  tropical  cyclones  by  tracking  clouds  using  successive  satellite  i  lanes. 
Their  results  are  encouraging  in  soite  of  some  difficulties  such  as  the 
short  lives  of  cloud  turrets,  subpixel  movement  of  clouds,  and  overcast  con¬ 
ditions  near  the  storm  centers.  An  experimental  oceanographic  satellite 
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known  as  SEASAT-1  during  its  short  lifetime  provided  an  additional  data 
source  to  define  large  and  mesoscale  wind  fields  near  tropical  cyclones. 

A  special  scatterometer  (SASS)  flown  with  SEASAT-1  measured  the  marine 
surface  microstructures,  and,  through  appropriate  algorithms,  marine  sur¬ 
face  winds  can  be  inferred. 

The  purpose  of  this  study  is  to  evaluate  the  impact  of  the  satellite- 
sensed  winds  on  tropical  cyclone  forecasts.  As  a  simulation  study,  data 
generated  by  numerical  models  will  be  used  in  place  of  real  data.  The 
general  strategy  of  such  simulation  studies  follows  that  of  Charney  et  al^ 
(1969).  First,  a  control  integration  of  the  numerical  forecast  model  is 
performed  to  generate  the  "true"  history  of  the  atmosphere  or  the  "observa 
tion".  A  series  of  "standard  forecasts"  is  then  generated  based  on 
different  initial  states.  Finally,  a  series  of  forecasts  with  the 
"observations"  assimilated  is  conducted  to  evaluate  the  impact  of  assimila 
tion.  For  detailed  reviews  of  such  simulation  studies  and  their  general 
strategies,  readers  are  referred  to  McPherson  (1975)  and  Bengtsson  (1975). 

As  a  preliminary  study,  an  axisymmetric  tropical  cyclone  model  is 
employed  in  this  study,  therefore  only  the  impact  on  intensity  can  be 
studied.  The  method  of  assimilation  used  is  the  dynamic  initialization 
by  relaxation  (DIR)  technique.  To  approximate  the  relationship  between 
the  real  atmosphere  and  forecast  models,  parameterized  physics  in  the 
model  that  generates  forecasts  (forecast  model)  are  altered  from  those 
in  the  model  that  generate  the  observations  (natural  model).  In  the 
following  sections,  the  numerical  model,  the  experimental  design,  and  the 
method  of  assimilation,  will  be  discussed  in  sequence.  Finally,  the 
results,  conclusions  and  proposed  future  research  will  be  presented. 
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2.  Numerical  Model 

The  axisymmetric  tropical  cyclone  model  used  in  this  study  is 
similar  to  the  one  described  in  Chang  (1977)  and  Anthes  and  Chana  (197G). 

The  governing  equations  are  in  primitive  form  and  are  in  o(=P/Ps)  coordinates. 
The  explicit  water  vapor  cycle  and  parameterization  of  cumulus  convection 
follows  Kuo  (1974)  and  Anthes  (1977).  The  boundary  layer  (BL)  is  contained 
in  the  lowest  model  layer,  parameterization  of  various  vertical  fluxes  is 
based  on  a  generalized  similarity  theory  in  which  Yamada's  (1975)  universal 
functions  are  used  (Chang  and  Madala,  1980).  Charnock's  equation  is  applied 
to  compute  marine  surface  roughness  length. 

The  model  atmosphere  is  divided  into  six  layers  (Table  1).  A  uni¬ 
form  horizontal  grid  interval  of  30  km  is  used  from  the  center  to  a  radius 
of  600  km.  The  grid  interval  is  progressively  increased  by  a  factor  of 
two  outside  600  km.  The  leapfrog  temporal  integration  method  with  the 
time-averaged  pressure  gradient  force  (Brown  and  Campana,  1978)  is  employed 
for  numerical  integration.  The  spatial  finite  differencing  is  of  the  second 
order.  The  mean  hurricane  season  sounding  (Sheets,  1969)  is  used  for  the 

initial  and  lateral  boundary  conditions.  The  coriolis  parameter  f  has  the 

-5  - 1 

constant  value  5  x  10  s  . 

A  period  of  36  h  of  the  control  run  during  which  the  model  tropical 
cyclone  underqoes  a  rapid  intensification  is  chosen  as  the  "truth"  or 
'observation"  (thereafter  referred  to  as  such  in  this  study  except  stated 
otherwise).  For  convenience.  -12  h  and  24  h  are  designated  as  the  start 
and  cei  of  this  oeriod.  Vie  select  a  nerioH  "f  rapid  intensification  for 
study  in  order  to  maanif.v  errors  in  the  forecasts. 
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Table  1 

Vertical  Structure  of  the  Model 


Layer 

| - 

t 

!  tPfrh)  for  Ps=  1000 

! 

Undisturbed  Height  (km) 
at  Layer  Center 

1 

0  -  200 

18.3 

2 

200  -  300 

10.6 

3 

300  -  600 

6.7 

f 

A  i 

600  -  800 

3.0 

5 

800  -  930 

1.2 

6 

930  -  1000 

0.3 

3.  Experimental  Design 

In  recent  years,  --any  simulation  studies  have  been  conducted  to 
evaluate  impact  of  incomplete  observation  data  on  numerical  predictions 
Charney  et  al,  1969:  Kasahara  and  Williamson,  1972;  Morel  and 
: 1  ?'Tsnd ,  1974;  Anthes,  1974;  Cane  et  aU  1S79).  In  a  similar  manner, 
numerical  integrations  conducted  for  this  study  can  be  grouped  into  three 
components  (Table  2): 

(1)  Nature  run  -  A  36  h  segment  of  life  history  of  tropical 
cyclones  designated  as  observation  as  defined  in  Section  2. 

(2)  Standard  forecasts  -  Two  36  h  forecast  starting  from  -  12  h 
an-  a  2^  h  forecast  starting  from  0  h.  The  initial  conditions 
for  standard  forecasts  are  obtained  by  the  static,  non- 
divergent  initialization  method  hased  on  the  nature  run. 
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(3)  Forecasts  with  assimilation  -  12  h  preforecast  integrations 
starting  fron  -  12  h,  during  which  satellite-sensed  winds 
are  assimilated  into  the  model  solution  followed  by  24  h 
forecasts  starting  from  0  h. 

A  unique  characteristic  of  previous  simulation  studies  is  that  a 
prediction  model  will  make  an  error-free  forecast  given  error-free  model¬ 
generated  initial  conditions  (Williamson,  1973).  Because  this  is  rather 
unrealistic,  errors  of  various  kinds  were  added  to  the  observations  either 
in  initial  conditions  for  the  forecasts,  or  in  data  for  assimilation. 

Both  random  errors  (e.g.,  Williamson  and  Kasahara,  1971)  and  bias  errors 
(e.g.,  Anthes,  1974)  have  been  introduced  into  the  observations  in  previous 
studies.  Forecast  runs  with  initial  random  errors  sometimes  exhibit 
unrealistic  error  growth  characteristics  because  gravity  waves  and  model 
physics  act  to  smooth  them.  Besides,  random  observational  errors  are  not 
the  major  problem  with  real  data,  where  systematic  errors  are  known  to 
have  caused  more  problems,  (McPherson,  1975).  Biased  errors  are  generally 
determined  subjectively  and  may  be  unwarranted  and  unrealistic.  No  error 
is  artifically  added  in  the  initial  fields  for  all  forecast  runs  in  our 
study,  instead,  errors  in  the  initial  wind  fields  are  introduced  by  the  non 
divergent,  gradient-balanced,  static  initialization  procedure  adopted  here. 
Such  initialization  procedure  is  currently  in  use  operationally.  Figure  1 
shows  the  errors  of  speed  in  the  initial  wind  fields  of  forecasts  at  -12 
h.  As  expected,  large  errors  occur  in  the  low-level  and  the  outflow  level, 
where  divergent  components  of  wind  vectors  are  largest.  The  initial  errors 
for  forecasts  initialized  at  0  h  are  the  same  characteristics. 
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List  of  Simulation  Experiments 


are  radial  and  tangential  winds,  respectively.  Subscripts  denote  layers. 


r(knt) 

Fig.  1  —  Wind  speed  errors  (m  s'-*-)  in  the  initial  <p  condition  at  -12  h  after  a  non -divergent,  static  initialization. 


In  previous  simulation  studies,  the  models  which  generated  the  fore¬ 
casts  were  identical  with  the  models  that  generated  the  observations  (see 
McPherson,  1975,  for  review).  This  of  course  is  very  unrealistic.  In 
reality,  numerical  forecasting  models  with  finite  spatial  resolutions  and 
parameterized  physics  cannot  reproduce  the  atmosphere  even  if  perfect  initial 
conditions  are  obtained.  To  properly  account  for  the  discrepancies  between 
real  forecasting  models  and  the  atmosphere,  the  parameterized  physics  in 
our  forecast  model  are  deliberately  altered.  The  parameters  changed  are 
those  we  feel  most  uncertain  about  in  the  physical  parameterizations  in  current 
numerical  models,  namely,  the  effective  air-sea  exchange  coefficients  and 
the  vertical  distribution  of  latent  heat.  The  effective  coefficients  of 
eddy  transfers  of  momentum,  sensible  heat  and  latent  heat  in  the  forecast 
model  are  set  at  90%  of  those  in  the  nature  model.  The  10%  error  is  well 
within  the  expected  error  in  the  BL  formulations.  The  vertical  distribu¬ 
tion  of  cumulus  heating  is  also  changed  so  that  approximately  5%  of  the 
heating  in  the  lower  troposphere  is  shifted  to  upper  troposphere.  The  5% 
difference  is  within  the  variation  of  the  observed  heating  distributions. 

Due  to  these  two  changes,  the  forecasts  in  our  study  do  not  asymptotically 
approach  the  observations  even  after  long  integration.  Because  the  errors  in 
the  parameterized  physics  in  our  forecast  model  are  within  the  differences 
between  the  atmosphere  and  the  current  operational  forecast  models,  the 
asymptotic  difference  between  our  "forecasts"  and  "observation"  is  quite 
realistic. 

4.  Method  of  Assimilation 

The  satellite-sensed  data  are  in  many  occasions  incomplete  in  that 
they  do  not  contain  observations  of  all  meteorological  variables  simul¬ 
taneously  or  the  observations  are  made  at  different  locations  and  times. 


To  incorporate  such  data  in  a  dynamically  consistent  way  into  the  numerical 
models,  suitable  methods  of  assimilation  must  be  used.  From  the  direct 
insertion  method  (e.g.,  Charney  et  a]_,  1969),  to  the  complicated  variational 
assimilation  (e.g.,  Sasaki,  1970),  there  are  many  methods  of  assimilation 
in  existence.  However,  not  all  methods  are  applicable  for  the  satellite- 
sensed  data  in  question.  The  wind  fields  derived  from  GOES  images  are 
basically  restricted  to  low  and  outflow-levels  in  tropical  cyclones 
(Rodgers  et  a[,  1977),  and  wind  fields  measured  from  SEASAT-1  are  at 
aneomometer  level.  For  such  data  with  poor  vertical  resolution,  a  method 
called  dynamic  initialization  by  relaxation  (DIR)  is  desirable. 

DIR  is  a  technique  wherein  the  meteorological  variables  are  relaxed 
(or  nudged)  by  using  the  model's  governing  equations  toward  the  observed 
values  during  a  preforecast  integration  (Anthes,  1974;  Hoke  and  Anthes, 

1976).  The  technique  has  shown  great  promise  in  real  data  applications 
(Nitta  and  Hovermale,  1967;  Davies  and  Turner,  1977;  Hoke  and  Anthes,  1977). 
Mathematically,  governing  equations  during  the  preforecast  integration  are 
modified  to: 

N 

-r  =  F  (X,t)  +  X  ^(t~n,  lSt>  >sr,  6z)  (x°  -  x)  (1) 

where  >'  is  an  element  in  the  vector  of  variables  X,  the  function  F  contains  the 
normal  terms  in  governing  equations,  x°  is  the  observation,  N  the  number 
of  observations,  and  \  the  relaxation  coefficient.  In  a  full  four  dimen¬ 
sional  assimilation,  A  is  the  function  of  observational  error,  cn,  the 
time  separation  of  the  observation,  5t,  the  horizontal,  (<5r),  and  the 
vertical,  (5z),  spatial  separations  between  observations  and  grid  points. 

It  should  also  depend  on  the  meteorological  variables. 


To  simplify  the  functional  form  of  A,  we  will  use  point-to-point 
relaxation,  i.e.,  variables  are  relaxed  toward  observations  made  at  the 
same  model  grid  points  only.  This  requires  that  observations  be  taken  at 
model  grid  points  and  all  observations  be  taken  simultaneously.  Note  that 
the  horizontal  resolution  of  the  satellite  measurements  do  approach  those  of 
typical  operational  forecast  model  of  tropical  cyclones.  With  the  develop¬ 
ment  of  suitable  BL  models  for  vertical  extrapolation  (Yu,  1980)  the  con¬ 
venience  of  point-to-point  relaxation  assumed  for  convenience  in  this 
study  is  nearly  available  in  operational  forecasting.  The  time  lag  o^ 
measurements  over  the  domain  of  tropical  cyclones  within  one  satellite 
revolution  is  negligibly  short  as  compared  to  the  12  h  "reforecast 
integration.  We  take  note  that  the  sv/ath  width  of  orbiting  satellites 
nevertheless  may  not  be  large  enough  to  cover  the  entire  tropical 
cyclone. 

The  satellite-sensed  winds  are  not  free  of  errors.  Rogers  et  al_ 

(1977)  estimated  the  mean  speed  errors  in  their  derived  winds  to  be 

2.5  m  s”1  relative  to  aircraft  measurements.  There  are  conflicting 

reports  on  the  errors  of  SEASAT  measurements  (Black,  1979;  Jones  and 

Piearson,  1978),  but  in  general,  the  errors  of  satellite-sensed  winds  are 

smaller  than  those  introduced  by  the  objective  analyses  over  the  oceans 

« 

(Cardone  e*  al_,  1976).  The  contribution  of  satellite  measurements  is  not 
in  the  general  error  reduction  but  in  the  filling  of  data-void  areas 
(Ghil  et  al_,  1979).  For  a  clear  demonstration  of  the  impact  in  assimilat¬ 
ing  winds  at  different  levels,  in  Exps.  2-9  it  is  justifiable  to  assume 
that  the  satellite  observations  are  error-free  in  comparison  to  the  initial 
and  model  errors.  However,  errors  of  different  magnitudes  are  added  to 
the  satellite-sensed  wind  in  Exps.  9E  and  9E2  to  evaluate  the  extent  to 
which  the  observation  errors  contaminate  the  forecast. 
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The  equations  of  motion  in  the  pre forecast  integration  in  this  study 
can  simply  be  written 

ft  *k  ■  F  <!•«>  *  »<Vk  -  »k)  (2) 

where  k  denotes  the  layer  in  the  model  where  observations  are  available. 

Three  different  values  for  A  are  tested:  Aw  =  10~4  s'1  for  weak  relaxation, 
-3  - 1 

Ac  =  10  s  for  strong  relaxation;  and  Aa  =  Ar  (6t  -  t)/12,  -12  h<  6t  0* 
for  attenuating  relaxation.  Figure  2  illustrates  the  time  variations  of  A. 
Observations  are  assumed  to  be  taken  at  0  h  and  are  assimilated  into  model 
prediction  during  -12  to  0  h  in  all  of  the  assimilation  experiments 
(Exps.  4  -  9). 


t(h) 


Fig.  2  —  The  three  relaxation  coefficients  used  in  DIRT:  X  w  for  weak 
relaxation,  Xg  for  strong  relaxation,  and  Xfl  for  attenuating  relaxation. 


5.  Standard  Forecasts 

Durina  the  oeriod  between  -12  h  and  24  h,  observation  shows  a  rapid 
intensification  of  the  trooical  cyclone,  the  minimum  central  pressure 
deepens  Trori  9r<8  mb  to  °53  mh  (Figure  3)  and  the  maximum  wind  speed  increases 
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from  29  m  s’*  to  52  m  s’*  (Figure  4).  A  36  h  forecast  starting  from  -12  h 
(Exp.  2)  and  a  24  h  forecast  starting  from  0  h  (Exp.  3)  are  conducted  based 
on  the  error-free  initial  mass  field  and  non-divergent,  gradient-balanced 
wind  fields.  Typical  wind  speed  errors  in  the  initial  wind  field  are 
illustrated  in  Figure  1. 

As  expected  from  previous  experience,  both  forecasts  have  an  initial 
dissipation  stage  due  to  the  onset  of  the  surface  friction.  The  weakening 
of  storm  intensity  is  especially  pronounced  in  Exp.  3  in  that  it  has 
larger  intensity  errors  before  12  h  than  Exp.  2  which  is  initialized  12  h 
earlier.  This  is  indicative  of  the  inadequacy  of  the  static,  non-divergent 
initialization  employed.  Improvement  of  forecast  during  the  initial  hours 
can  be  achieved  by  using  a  divergent  static  (Tarbell,  1979)  or  a  dynamical 
(Hovermale  and  Livezey^W?;  Kurihara  and  Bender,  1979)  initialization 
scheme. 

After  the  radial  circulations  develop,  both  forecasts  reproduce  the 
observed  intensification  but  at  slower  rates.  After  12  h,  the  24  h  fore¬ 
cast  (Exp.  3)  yields  better  prediction  than  the  36  h  forecast  by  approximately 
2  mb  in  minimum  pressure  and  2  m  s’*  in  maximum  wind  speed.  Both  forecasts 
predict  weaker  storm  intensity  as  compared  to  the  observation.  The 
difference  between  the  observed  and  predicted  intensities  at  24  h  is  about 
15  mb  in  minimum  pressure  and  10  m  s’*  in  maximum  wind  speed.  The  divergence 
of  the  forecasts  from  the  observation  is  a  consequence  of  the  "imperfect" 
physical  parameterizations  in  the  forecast  model. 

We  select  the  root-mean-square  errors  (e)  as  a  measurement  of  the 
accuracy  of  the  predictions  (Panofsky  and  Brier,  1968).  Evolutions  of  5 
for  wind  speed  (V),  temperature  (T),  and  specific  hunidity  (q)  with  respect 
to  the  observation  for  Exps.  2  and  3  are  shown  in  Figures  5,  6,  and  7, 
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respectively.  The  initial  e(V)  is  large  at  about  4  m  s"^  in  both  experi¬ 
ments  due  to  the  non-divergent  initialization.  It  decreases  for  the  first 
18  h  in  both  experiments  as  the  forecast  storms  intensify.  As  evident  in 
Exp.  2  after  6  h,  values  of  e( V)  begins  to  increase,  indicating  a  deteriorat¬ 
ing  forecast. 

e(T)  and  e(q)  increase  rapidly  after  initialization  with  time  from 
the  error- free  mass  field.  Their  values  escalate  to  1.5°K  in  temperature 
and  0. 9  —  1.1  g  kg  ^  in  specific  humidity  at  24  h.  As  expected,  Exp.  3 
produces  a  better  prediction  than  Exp.  2  during  most  of  the  period  0  -  24  h. 

6.  Forecasts  with  Assimilation  of  Low-level  Winds 

In  Exp.  4,5,  and  6,  the  observed  low-level  radial  (ug)  and  tangential 
(Vg)  winds  at  0  h  are  assimilated  by  DIRT  into  the  24  h  forecast  during  a 
pre-forecast  integration  from  -12  to  0  h  (Table  2).  The  relaxation  coef¬ 
ficients  are  \W,XS,  and  Xg  in  Exps.  4,  5,  and  6,  respectively  (Figure  2). 

Figures  8  and  9  show  the  minimum  pressures  and  the  maximum  wind 
speeds  for  these  three  forecasts  with  low-level  winds  assimilated.  It  is 
apoarent  that  the  DIRT  with  the  weak  relaxation  coefficient  (Exo.  4)  does 
not  alter  the  prediction  appreciably  toward  the  observation  durinq  the  pre¬ 
forecast  integration.  The  followinq  24  h  forecast  has  no  apparent  improve¬ 
ment  over  the  standard  forecasts. 

The  model  adjustments  are  considerable  when  strong  and  attenuating  relaxa¬ 
tion  coefficients  are  applied  in  Exps.  5  and  6.  The  maximum  wind  speed  in 
the  oreforecast  integration  converges  to  the  observed  value  within  a  couple 
of  hours.  The  minimum  pressure  also  approaches  the  observed  value  at  0  h 
within  hi  h,  in  agreement  with  the  theory  of  geostrophic  adjustment.  The 
maximum  wind  soeed  in  Exp.  5  achieves  the  observed  value  due  to  constantly 
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strong  relaxation.  However,  as  the  assimilation  terminates  at  0  h,  strong 
model  adjustments  occur  in  both  experiments.  The  model  rejection  occurs  in 
Exp.  5  where  the  maximum  wind  speed  decreases  approximately  5ms  *  in 
three  hours  and  stays  lower  than  that  of  the  standard  forecast  Exp.  2.  The 
rejection  is  similar  in  Exp.  6.  so  there  is  still  no  improvement  in 
intensity  forecast. 

we  thus  conclude  that  low-level  wind  observations  are  not  beneficial 
to  intensity  forecasts  of  tropical  cyclones  if  assimilated  by  DIRT.  We 
can  also  conclude  that  the  attenuating  relaxation  coefficient  is  more 
effective  in  assimilating  the  observed  data  (cf.  Fig.  9)  and 
desirable  for  eliminating  model  adjustments. 

However,  an  examination  of  forecast  errors  in  Exps.  4,  5, 
and  6  is  warranted.  During  the  preforecast  integration,  e(V)  decreases 
% -i tli  time  as  relaxation  forces  the  low-level  wind  to  asympotically  approach  the 
0  h  observation.  The  e(V)  value  from  Exp.  5  at  0  h  reaches  the  lowest  level 
of  all  (Figure  10).  The  model  adjustments  cause  the  error  to  be  at  levels 
higher  than  those  of  the  standard  forecasts  (Exp.  2  and  3)  after  6  h  since 
the  assimilation  has  been  rejected  by  the  model.  The  e(T)  and  e(q)  are 
similar  to  e( V)  in  that  they  decrease  with  time  in  the  preforecast  integration 
for  Strong  relaxation  and  they  subsequently  increase  to  levels  equivalent 
or  higher  than  those  of  the  standard  forecasts. 

The  rejection  of  assimilation  of  low-level  winds  in  above  experiments 
can  be  attributed  to  the  insufficient  vertical  coupling  between  the  lo'-- 
level  an"'  high-level  momentum  fields  during  the  12  h  period  of  Dreforecast 
integration.  A  longer  period  of  preforecast  integration  may  produce  enough 
vertical  coupling  thro.'1-!-1  ■"odel  dynamics  and  oh.vsics,  but  is  no*  very 
meaningful  in  practice.  It  is  then  logical  to  test  assimilatior  of  additional 
winH  observations  at  higher  levels  since  they  can  be  made  available 

(Rodger  et  a]_,  1977) . 


maximum  wind  speeds  of  the  observation  (Exp.  1)  and  standard  forecasts 


t(h) 

Fig.  7  —  The  time  series  of  the  forecast  rms  errors  in  specific  humidity  for  Exps.  2  and  3. 


Fig.  9  —  The  maximum  wind  speeds  of 


7. 


Forecasts  with  Assimilation  of  Low  and  Higher-level  Winds 

As  listed  in  Table  2,  observations  of  higher  level  winds  in  addition 


to  t'-''  low-level  winds  are  assimilated  by  DIRT  using  relaxation  coe¬ 
fficient  \  .  In  Exp.  7,  v r  (check  Table  1  for  pressure  level)  at  0  h  are 

cl  t) 

assimilated;  in  Exp.  8,  v^  and  v^;  and  in  Exp.  9,  U2  and  v^. 

The  assimilations  of  higher  level  winds  yields  significantly  different 
intensity  forecasts  from  the  standard  forecasts  as  evident  in  Figures  11 
and  12  showing  the  minimum  pressure  and  the  maximum  wind  speed,  resoectively . 
In  addition,  the  three  experiments  forecast  very  different  minimum  pressures 
even  during  the  preforecast  integration  where  the  same  X  is  used.  Amono  the 
three,  Exp.  8  yields  the  best  prediction,  with  maximum  difference  of  only  4.5 
mb  in  central  pressure,  at  24  h.  Exp.  7  predicts  a  very  intensive  storm,  with 
central  pressure  deepening  to  955  mb  at  12  h.  Exp.  9  predicts  a  weaker 
storm  than  the  observation,  however,  it  forecasts  better  than  the  standard 
forecasts  (Exp.  2  and  3)  and  forecasts  where  only  low-level  winds  are 
assimilated  (Exps.  4,  5,  and  6). 

It  is  interesting  that  the  storm  intensities  in  Exps.  7  and  8  are 
drastically  different  when  the  only  difference  in  the  experiments  is  that 
the  observations  of  v^  at  0  h  are  available  for  assimilation  in  Exp.  7.  As 
demonstrated  by  the  vertical  profile  of  v  at  r  =  30  km  in  Figure  13,  there 
is  a  strong  vertical  shear  in  the  0  h  tangential  wind  observation.  Tangen¬ 
tial  wind  speed  decreases  upward  associated  with  the  strong  warm  core  at 
level  4.  In  Exp.  7  the  too  strong  storm  intensity  is  due  to  the  assimilation 
of  the  stronger  circulation  below  800  mb,  whereas  the  vertically  decreasing 
tangential  circulation  and  effects  of  warm  core  are  properly  assimilated  in 
Exp.  8.  This  suggests  that  when  observational  data  are  to  be  vertically 
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interpolated  in  diagnosis  or  analysis,  strong  vertical  shear  and  the 
related  baroclinic  effect  must  be  taken  into  account. 

The  forecast  errors  in  wind  speed,  temperature  and  water  vapor  are 
shown  in  Figures  14,  15,  and  16,  respectively.  The  values  of  i(V)  decrease 
in  time  during  the  pre forecase  integration,  as  in  Exps.  4,  5,  and  6.  Note 
the  e(V)  in  Exp.  9  is  the  smallest  because  initial  errors  are  largest  at 
the  assimilated  levels  (Figure  1).  The  errors  in  Exps.  8  and  9  remain 
smaller  than  those  of  the  standard  forecasts,  especially  in  Exp.  9,  where 
the  error  is  50%  lower. 

The  error  in  temperature  field  of  Exp.  7  arises  early  in  the  prefore¬ 
cast  integration.  This  shows  that  the  effects  of  the  warm  core  on  tangential 
circulation  are  not  properly  assimilated  as  mentioned  earlier.  The  errors 
in  Exps.  8  and  9  are  generally  smaller  than  for  the  standard  forecasts 
throughout  the  24  h  forecast  period.  The  specific  humidity  errors  for  these 
three  forecast  experiments  are  higher  than  for  the  standard  forecasts  with 
Exp.  7  having  the  highest  error.  Between  -12  and  6  h,  e(q)  in  Exp.  8  is  very 
low  because  the  inflow  and  outflow,  which  nearly  determine  the  net  total 
water  vapor  convergence,  are  assimilated. 

The  higher  e(q)  in  Exps.  8  and  9  in  spite  of  the  better  intensity 
forecasts  and  lower  i(V)  and  e(T)  is  probably  due  to  the  different  physical 
parameterizations  in  the  forecast  model. 


Fig.  11  —  The  minimum  pressures  of  the  observation  (Exp.  1),  and  standard  forecasts  (Exps.  7,  8,  and  9). 


12  -6  0  6  12  18  24 
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.  16  —  The  time  series  of  the  forecast  rms  errors  in  specific  humidity  for  Exps.  7,  8,  and  9. 
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8.  Deterioration  of  Forecast  due  to  Satellite  Observation  Errors 

In  previous  sections,  we  have  examined  the  impact  of  assimilating  the 
error-free  satellite  winds.  The  assumption  of  perfect  satellite  observation 
is  for  clearity  in  comparison.  Satellite-sensed  winds  are  of  course  not 
error-free.  As  discussed  earlier,  the  mean  speed  error  in  satellite  winds 
vary  from  2.5  m  s'1  (Rodgers  et  aj_,  1977)  to  as  large  as  8  m  s'1  (Black, 
1979).  We  now  turn  our  attention  to  the  influence  of  the  satellite  observa¬ 
tion  errors  on  the  forecast. 

Because  low-  and  outflow  level  winds  are  most  likely  to  be  obtainable 
operationally,  we  repeat  Exp.  9  with  artificially  introduced  observation 
errors  in  both  the  initial  condition  and  satellite-sensed  winds.  Randomized 
error  of  about  2.5  m  s'1  and  biased  errors  are  added  into  the  observation 
at  -12  h.  The  biased  errors  have  a  maximum  of  8  m  s'1  at  r  =  30  km  and 
decreasing  with  radius  to  zero  at  r  =  150  km.  These  errors,  after  the 
balanced,  static  initializations  are  equivalent  to  approximately  0.5°K 
random  errors  in  temperature  field  and  4  mb  error  in  central  pressure.  Four 
experiments,  Exps.  2E,  3E,  9E,  and  9E2  are  carried  out  based  on  initial 
conditions  containing  such  errors.  The  satellite-sensed  winds  at  0  h  con¬ 
tain  random  errors  with  maximum  speed  error  of  2.5  m  s  '  in  Exp.  9E  and 
5  m  s'1  in  Exp.  9E2.  Exps.  2E  and  3E  are  identical  to  Exps.  2  and  3  except 
for  the  introduced  errors  in  the  initial  condition. 

As  summarized  in  Table  3,  the  average  12  and  24  h  forecasts  in  Exps. 

2E,  3E,  and  9E  are  worse  than  their  respective  counterparts  in  error-free 
simulations.  For  example,  Exp.  9  has  an  averaged  forecast  error  of  11  mb 
in  minimum  pressure,  6.5  m  s'1  in  maximum  wind  speed,  and  2.5  m  s'1  in  e(V), 
whereas  Exp.  9E  has  an  averaged  forecast  error  of  16  mb,  11.7  m  s’1,  and 
4.5  m  s'1.  It  is  encouraging  that  Exp.  9E,  in  which  the  magnitude  of  the 
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observation  errors  are  typical  for  operational  forecast,  is  a  better 
forecast  than  the  standard  forecasts  of  Exps.  2E  and  3E.  It  forecasts 
better  than  the  standard  forecasts  by  4  mb  in  minimum  pressure,  approxi¬ 
mately  2  m  s”1  in  maximum  wind,  and  0.5  in  i(V).  It  is  also  interesting 
that  Exp.  9E2  performs  only  slightly  worse  than  Exp.  9E,  although  the 
error  level  i:>  twice  as  large  as  that  in  Exp.  9E. 

From  the  above  comparison,  we  conclude  that  the  errors  in  satellite 
winds  could  lead  to  a  deteriorate  forecast,  and  that  the  assimilation  of 
satellite  low-  and  outflow-level  winds  can  improve  the  forecast  if  these 
errors  are  less  or  equal  to  those  contained  in  the  initial  wind  field. 
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9. 


Summary  and  Discussion 

The  impact  of  accurately  measured  marine  surface  winds  of  sufficient 


spatial  coverane  and  resolution  on  the  24  h  intensity  forecast  of  tropical 
cyclones  'een  studied  with  simulation  experiments.  The  model  physics 
in  the  forecast  model  were  altered  from  those  in  the  nature  model.  The 
observations  are  assimilated  into  the  numerical  forecast  with  dynamical 
initialization  by  relaxation  during  the  pre-forecast  integrations  from 
error-free  mass  fields. 

The  results  indicate  no  improvement  in  forecast  accuracy  when 
low-level  winds  are  assimilated  according  to  the  abovementioned  pro¬ 
cedure.  We  note  that  a  strong  relaxation  coefficient  causes  rejection 
of  the  assimilation  within  a  few  hours  of  forecasting  and  that  a  weak 
relaxation  coefficient  is  ineffective. 

Significant  improvements  are  achieved  when  all  winds  below  600  mb 
re  assimilated.  This  conclusion  is  easily  understandable  because  the  major 
characteristics  of  the  tropical  cyclone  such  as  the  vortex  strength,  the 
warm  core,  and  the  vertical  shear  are  included  in  such  observations.  But 
simultaneous,  high-resolution  observation  required  for  such  assimilation 
is  very  difficult  to  obtain.  It  is  encouraging  that  improvement  in  forecast 
can  also  be  achieved  when  low  and  outflow-level  winds  are  assimilated 
because  wind  fields  at  these  two  levels  are  most  likely  available  from 
satellite  observations. 

The  forecast  with  low  and  outflow-level  winds  assimilated  worsens  with 
increasing  observation  errors.  However,  even  if  the  root -mean-square 
error  in  the  satellite  observation  is  equivalent  to  that  in  the  initial 
wind  field,  assimilation  of  low  and  outflow-level  winds  still  improves  the 
forecast . 
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Caution  must  be  taken  in  interpreting  these  findings  for  operational 
applications,  as  is  the  case  for  all  simulation  studies,  because  the  extent 
to  which  they  approximate  reality  is  difficult  to  determine.  The  finding 
that  the  low-level  wind  observations  alone  cannot  improve  the  forecast  when 
assimilated  by  DIR  should  not  cast  doubt  on  the  usefulness  of  observing 
systems  which  measure  marine  surface  winds.  Since  an  axisymmetric  tropical 
cyclone  is  employed  in  this  study,  the  position  of  the  storm  is  assumed 
known.  Also,  the  mass  fields  are  assumed  to  be  error-free  in  the  static 
initializations  of  the  forecasts.  The  precise  center  location  and  perfect 
mass  field  are  not  commonly  available  for  operational  forecasts,  where 
meteorologists  have  to  be  content  with  uncertainties  of  the  storm  center 
position  and  with  the  "bogussed"  circulations.  The  marine  surface  winds 
are  invaluable  in  defining  the  low-level  circulations  and  in  locating  the 
storm  centers  which  otherwise  would  be  impossible  over  data-void  oceans. 

Since  our  results  with  assimilation  of  the  low  and  outflow-winds  are 
encouraging,  and  since  these  will  be  the  focal  levels  in  satellite  observa¬ 
tions,  future  research  with  a  three-dimensional  tropical  cyclone  model  is 
warranted.  In  a  three-dimensional  study,  the  impact  of  the  satellite- 
sensed  winds  on  storm  track  forecast  can  be  investigated.  The  effects  of 
time-lag  within  one  satellite  revolution  discussed  earlier  and  the  effects 
of  the  swath  width  can  also  be  studied.  Finally,  real  data  case  studies 
can  be  carried  out  with  a  three-dimensional  model. 
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